Theory of multiwave mixing and decoherence control in qubit array system 
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We develop a theory to analyze the decoherence effect in a charged qubit array system with photon 
echo signals in the multiwave mixing configuration. We present how the decoherence suppression 
effect by the bang-bang control with the n pulses can be demonstrated in laboratory by using a bulk 
ensemble of exciton qubits and optical pulses whose pulse area is even smaller than tt. Analysis 
is made on the time-integated multiwave mixing signals diffracted into certain phase matching 
directions from a bulk ensemble. Depending on the pulse interval conditions, the cross over from 
the decoherence acceleration regime to the decoherence suppression regime, which is a peculiar 
feature of the coherent interaction between a qubit and the reservoir bosons, may be observed in the 
time-integated multiwave mixing signals in the realistic case including inhomogeneous broadening 
^ . effect. Our analysis will successfully be applied to precise estimation of the reservoir parameters 

O ' from experimental data of the direction resolved signal intensities obtained in the multiwave mixing 

I technique. 

CN ■ PACS numbers: 03.67.Pp, 42.50.Md, 78.20.Bh 



■ I. INTRODUCTION 

o 

1—1 ■ 

The analysis of decoherence of elementary excitations in solid state is one of the main approaches to understand 
physics in condensed matter. It also attracts much attention in device physics attempting to implement quantum 
information processing devices. A central issue is how to realize and control coherent evolution of well coded quantum 
(-H states. Effective means should be chosen depending on the way of coupling to the system and decoherence mechanism. 
^ I' Among various proposals semiconductor devices have potential advantages, including the ease of integration and the 
^ ' , , use of matured industrial technologies. 

P5 ' In particular, excitons in semiconductor make clear distinction from other two level systems for qubits in the sense 
, that it can efficiently couple to photons, which are the signal carriers to build a communication network. In addition, 
^ its large dipole moment makes possible direct optical control of qubits in a time scale of femto second. Various results 
'[ of optical coherent manipulation of exciton in semiconductor have been reported, such as indirect P, |2| and direct 
^ ■ [3 observations of excitonic Rabi oscillations in higher dimensional semiconductor structures, manipulations of a 
i one-qubit rotation of a single quantum dot exciton 0, S IE S > direct observations of excitonic Rabi oscillations in 
i^jj quantum dot ensemble and entanglement manipulation in quantum dots jlCtll ij . They open great opportunities 
for quantum information processing based on exciton qubits. 

On the other hand, exciton qubits always suffer from unwanted coupling with the external degrees of freedom 
of charged excitations including the phonon scattering and the scattering by excitons themselves, leading to the 
decoherence. In particular, the fast decoherence time of excitons which is typically a few tens of pico second will be 
the first obstacle to be overcome. 

The well known methods to fight the decoherence are quantum error correcting codes [T^ IT^ [T3 . Iisl IT^ IT^ and 



encoding with decoherence free subspace [18|, ll9|, |20|, |21|, |23|, |2J| . Both methods rely on encoding with multipartite 
entangled states for correcting errors by redundant qubits, or for making the state evolution immune to decoherence 
by exploiting symmetries of the system-environment interaction. It is yet very challenging to prepare the multipartite 
entangled states in the desired form in solid state device. 

The other, and only plausible method with present technologies might be the dynamical decoupling control by 
driving the system continuously with external pulsed field sequences, the so called bang-bang (BB) control 25, 2a.l27l 
|2l|29|,|3d,|3ll|32|,(3l|MIM|3il^, or with p ro j ecting onto a certain measurement basis (38|, |39|, |40j , the so called 
quantum Zeno effect (for recent review, see |4lll42l |'). This does not basically require complicated qubit encoding, but 
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effective at a single qubit level. The former is easier from the technical viewpoint, and therefore to be demonstrated 
as the first step for solid state qubits. The BB control relies on a time reversal of the decoherece process in a short 
time scale within the reservoir correlation time. Similar ideas have been known such as spin echoes |43| , photon echoes 
[4^ , wave packet echoes 'iE] , and charge echoes 0| . These techniques have been powerful tools to investigate 
transient phenomena. The BB control can be regarded as a variant of them, focusing on an active decoherence control 
of qubit. This has been demonstrated recently with a solid-state nulear qubit [48l |. 

Our concern here is the BB control with optical pulses for exciton qubits in semiconductor. In order to put this 
method to the experimental test, we will face mainly two problems. One is the difficulty of detecting the expected 
stabilization effect by weak signals from a single exciton qubit. The other is that the presently available laser power is 
not strong enough yet for preparing the tt pulse area. One possible solution is to use a bulk ensemble of excitons. Strong 
signals from the ensemble make the detection easy. Furthermore they make it possible to observe the decoherence 
suppression effect due to the ideal tt pulses even if the irradiated pulses have the area of smaller than tt. This is 
because some fraction of qubits can effectively feel the ideal tt pulses, and they form a macroscopic polarization lattice 
which radiates photon echo signals in a certain phase matching direction. The bulk ensemble is, however, always 
associated with irrelevant phenomena to the present issue, such as inhomogeneous broadening of the energy levels 
and complicated light diffraction due to polarization lattices of excitons. 

It is then necessary to extend the theories for a single qubit to a qubit array system to analyze the effect of our 
interest with rather complicated signals. The purpose of this paper is to develop a theory to analyze the decoherence 
with the photon echo signals in the multiwave mixing configuration, and to predict the effects of decoherence control 
in a practical experimental setting. 

The paper is organized as follows. In section^ we present the model and introduce basic notations. In section UTTI 
we first develop a formulation to evaluate the multiwave mixing signals in the ideal tt pulse case. We then present 
the results in the case of homogeneous broadening, and introduce basic notions to analyze the state evolution of the 
qubit-reservoir system using the multiwave mixing signals. In section Hvl we develop a formulation for the weak pulse 
case, and present numerical simulations on the time-resolved and time-integated multiwave mixing signals. They 
include not only the suppression but also the acceleration of decoherence induced by the optical BB control. Section 
Ivlis for concluding remark. 

II. MODEL 

In this paper we consider an ensemble of qubits and identical reservoirs. The model Hamiltonian is given by 

H{t)=Ho+HQR + HQFit), (1) 



where 



i;iz.,<7,(j) + E?ir!,st(j)B,(j) , (2) 

Hqr = E (^■) E [ai^l (■?■) + 9tBi ij)] , (3) 
H^At) = E £^m(Oe'"'"""^^"*l + c.c. (4) 

j m 

Hq describes the free evolution of the qubit and the reservoir systems. The index j refers to the jth site of qubit 
whose ground and excited states are | l)j and | t)j, respectively. Uj is the energy separation of the jth qubit. CTz is 
the Pauli z component operator. We assume that the qubit at the jth site interacts with its own reservoir around 
the site, and that the reservoirs at different sites are independent but have an identical mode structure {^1} and a 
temperature T. The index I is the Ith mode of the reservoir boson with the angular frequency fi; which is described 
by the creation and annihilation operators B\{j) and Bi{j). Hqr describes the coupling between the qubit and the 

reservoir systems by the coupling constant gi. HQF{t) describes the qubit-external field interaction in the rotating 
wave approximation, k is its coupling constant. CT-|_(j) acts as (T_|_(j)| J,)^ = | Vj is the coordinate of the jth site 
of qubit. Em{t), k^, and uj are the external (classical) field amplitude, its wave vector, and angular frequency of the 
mth pulse. 
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FIG. 1: Optical pulse sequence. 



The essential physics of the decoherence in the above model is the same as that in the model consis ting of a single 
qubit and a bosonic thermal reservoir, and has already been clarified by extensive theoretical analyses |l8ll25ll26ll37l| . 
The new feature included in the above model is just an array structure of qubits. Main concern here is to know 
how the decoherence occuring at a single qubit level can be detected in the multiwave mixing signals from this array 
structure. 

The interaction pictured Hamiltonians are given by 



j I 



(5) 



(6) 



For simplicity, we assume that the optical field in Hqp (t) are composed of a sequence of optical pulses with rectangular 
temporal shape as shown in Fig. ^ which are specified by 



Emit) 



otherwise 



(7) 



We further assume that the pulse duration r™ is so short that the effect of the qubit-reservoir interaction HQ^iit) can 
be neglected during the period [tm,tm + T,n]- Taking the limit of — > and Em oo with the finite pulse area 
6m = —'^iEmTmi^, the State evolution under the mth optical pulse is described by the unitary operator 



(m) 
QF 



n 



~(m) 



(8) 



with 



~(m) 



where '/'j'"'' ~ k,„ • Tj + (vj — Lo)tm- Each qubit evolves according to 

i>J = cos^l i>j +e'"^^- 'sin^l t)i, 



(9) 



(10) 



^fffi l T)j = cos^l ])j - e 'sin^l 
The state evolution between the pulses are described by the qubit-reservoir interaction 



(11) 



(m) 



QR 



Texp 



dt'HQnit') 



n 



(12) 
(13) 
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with 



~(m) 



where 



and 



0(im+li*ni) — E 1^' 



(m) _ 9i r _ „jn,t™+i] 



To evahiate the state evolution by Uq^ , we may use the formula 

-(m) 



«SlT),®l*>R,,-|T>,®^5^l*>R„ 

where the state |\E')Rj represents any pure state component of the jth reservoir and 



with the displacement operators 



i^,(±aj™)) = exp (±ari?/(j) T ar*S/(j) 



(14) 

(15) 
(16) 

(17) 
(18) 

(19) 
(20) 



For the state |5')rj, we will later substitute the Fock state basis \ni...ni...)j which describes the occupation of the 
reservoir boson modes. 

The total evolution operator for the M pulse sequence is given by 



(21) 



with 



As for Uq^j, we use 



~(M)^~(M} ^~(M-1)~{M~1) ~(0) -.(0) 



QRj "QFj "QRi "QFj 



"QRj'"QFj"QRi'"QFj- 



(22) 



(23) 



i.e., is understood as t, instead of the definition of Eq. llt)|) . The final state of the total system is represented 

by 



where we assume the initial state like 



pQR(io)=nu).^(ii®n/^R^v(^) 



(24) 
(25) 



with the thermal state 



E^ (i+^',)«,o-)+i '"'(^')^^^-<"'(^')'' 



n, (i)=0 



(26) 
(27) 
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where 



ni = 1/ (^e'^^'/'^BT _ 
By using the formula Eqs. p7|l and (|18|l . the final state can be represented as 



(28) 



PQR 



where 



\PnU){t))QRj = Uj{t,to)\ l)j (g) \n{j))nj, 



(29) 



(30) 



and n(j) = (no ( j ),..., n;(j), ...) denoting the boson number occupation. To describe the state of Eq. H3()|l . we introduce 



the two operators Fj^^ and G^-*^^ by 



|p^(,.)(t))Q«, ^ (-1)* ( I i), ® Gf ) + I T), ® ^^^"^ ) |n(j-))R., 



for M =even, and 



for M =odd. 



|P«0) W)Qi?j = (-!)" 



-U),®F,^'"^ + |T),®Gf))|n(j))R„ 



(31) 



(32) 



III. STATE EVOLUTION UNDER n PULSES 



A. Formulation 



In this section we consider the pulse sequence consisting of the 0th exciting pulse with the area = '^/'^ a-t time 
to, and the subsequent M 7r-pulses, i.e., 6m — tt (m — 1, ...,M). By using Eqs. (|10|l . (|17|l . and (|18|l . we obtain 



p(A/) _ J_ »(<[.f^'+0f>)f,(M) 



and 



where 



G 



(M) _ 1 A(Af) 



V2 



(33) 



(34) 



M 



(35) 



The decoherence property of the qubit system is measured by observing the intensity of the radiation from the 
macroscopic polarization. The polarization operators in the interaction picture are defined by 



5±(q,<)^^a±(j)e^'i-"-^±''^- 



(36) 



The macroscopic polarization is then evaluated by 



P(q,t) = Tr 



^5+(q, t)pQ^{t) 
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, (M = even), 



, (A/ = odd), 



(-1) 



M 

-E 



expi 



q • r 



xTrr 



. I \ m=0 / 



Let us introduce the decoherence exponent by 



-r(t) 



Tr 



R 



M 



The decoherence exponent, which is now independent of j, is expressed as |18l \2t 

2 



r 



M 



■m=0 



coth 



= 2^^/(f7,,i)coth 



2fcBT 



where 



and 



M-l 



By substituting the definition of (Eq. (123) and 0^™'' = k„j • Vj + [vj — uj)tra into Eq. H37() . we have 



^^(q,t) 



(-1) 



M 



-r(t) 



^expi 



(q - (-l)^'^K(^^)) . r, + (z., - u;) Am + (-1)^^ ^ (-1)™A„ 



m=0 



where 



M 



K(^^) =2^(-l)"k„+ko, 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



Am = tm+i — tm, and Am — t — tm- Eq. (|42|l shows that the radiation signal is observed in the phase matching 
direction q= (-1)*^K(*^). 



B. Homogeneous case 



We first consider the homogeneous case of = ^^.^for aU i. This case was already studied by Aihara as the photon 
echo in the four wave mixing (two pulse) scheme IH^ . He pointed out that the interaction between the qubit and 
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FIG. 2: The time-resolved multiwave mixing signal intensities under the multi n pulse irradiation (solid lines) compared with 
the ones of the free induction decay (one dotted lines) in the case of the ohmic reservoir model. The parameters are f2c=8meV, 
a = 0.1 and T^IOK. M is the number of pulses and A is the pulse interval. 
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FIG. 3: The time-resolved multiwave mixing signal intensities under the multi n pulse irradiation (solid lines) compared with 
the ones of the free induction decay (one dotted lines) in the case of the Gaussian-ohmic reservoir model. The parameters are 
r2p = 13meV, 7p=4meV, Op — 0.05, f7c=8meV, a = 0.1, and r=10K. M is the number of pulses and A is the pulse interval. 

the reservoir bosons induces the slow frequency modulation on the qubit in the short time region where the reservoir 
keeps the phase memory, and plays the similar role to the inhomogeneous broadening. In fact he predicted that the 
photon echo must be observed even in the homogeneous case. The BB control proposed by Ban [23, and Viola and 
Lloyd 26] can be regarded as a natural extension of this principle to the multi-pulse scheme. It is well studied how 
to design a pulse sequence for effective suppression of decoherence I^B, 113 ■ The decoherence can be made worse 
by applying a pulse sequence, if the pulse intervals are not taken short enough. 

Fig. [5]shows numerical results of the decoherence under an M tt pulse sequence with an interval A. The solid lines 
represent the time-resolved multiwave mixing signal intensities of the multiwave mixing signals, Iq{t) = |P(q, i)p, in 
the phase matching direction q = (— i)*^k(*^) (Eq. (|43|l ). while the one dotted lines represent the ones of the free 
induction decay signals in the direction kg. We have assumed the ohmic reservoir model characterized by the spectral 
density 



= anexp{~ — ), 



(44) 
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where fie is the cutoff frequency and a is the dimensionless couphng constant. As seen, for a pulse interval of 
A = 0.2 ps, the decoherece is not suppressed but accelerated. In Fig. O we have shown the same kind of results in 
the case of the Gaussian-ohmic reservoir model, which has not only the ohmic character with a wider spectrum but 
also the Gaussian character with a sharper spectrum around a characteristic frequency Qp, such as a specific phonon 
mode. The decay profile exhibits the oscillation characterized by the frequency Qp. For shorter pulse intervals, the 
decoherence is accelerated by adding tt pulses. 



I{n) = afiexp(-— ) 



exp 



7p^ 



(45) 



Although this decoherence acceleration phenomenon is essential to understand the decoherence, its interpretation 
given so far is unsatisfactory. Here we revisit this point by considering a simplified model 



HQKM)=ha,igpB^ +g;B), 



(46) 



i.e., a signal qubit interacts with a single mode field of the reservoir boson oscillator with the angular frequency fip. 
Under this Hamiltonian, the qubit couples with the reservoir as 



(I T> + U))®|n)R 

t> <E) D[a{t - to)] + \i)<E} D[-a(t - to)] ) |n)R, 



where |n)Fi, is the n boson state component, and 

i) {±a{t - to)) = exp ±a{t - to)B'' T a(t - to)* B 

with 

a(t-to) = ^(e'f^p*o -e'"''*). 
The free induction decay is characterized by the decoherence exponent 

r.(t)^2M;^(,)eoth(|^).,'/.(0, 



(47) 



(48) 



(49) 



(50) 



where 



= 2< 1 - cos 



Vlpit-to) 



(51) 



Thus r+(t) oscillates in time and the qubit is recohered, i.e., disentangled from the reservoir, at t — io — 
(27r/r2p)x (integer). On the other hand, when a vr pulse is applied at a time ti, the decoherence exponent is given by 



with 



r_(t)^g'/-(0 



(52) 



= 2[3-2cos(fipAo)] 

— 2y^5 — 4cos(r2pAo)cos[fip(t — to) + ^p], 



where 



For small enough Aq (= ti — to), 



t&mp 



2cosf2pAo 
1 — 2cosfipAo 



(53) 



(54) 



r_(t) ~2g'{l - cos [fip(t-to-2Ao)]}. 



(55) 
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FIG. 4: The time-resolved multiwave mixing signal intensities under the multi n pulse irradiation (solid lines) compared with 
the ones of the single tt pulse, i.e., M = 1 (one dotted lines) in the case of the ohmic reservoir model when the inhomogeneous 
broadening effect is taken into account by setting 5B=5meV. The other parameters are the same as in Fig. |2| 



This should be compared with r+(i). We can see that the dephasing exponent r_(f) is the one shifted in time by 
2Ao from T+{t). Namely, by applying the tt pulse soon after the state | t) + | 4) is prepared, the decoherence starts to 
take place at the time 2Aq later compared to the free induction decay without any tt pulse. For larger Ag, however, 
this is not true. As an extreme case, we take fipAo ~ tt. Then 

T-W = g' {6cos [rjp(t - to)] + 10} > r+{t). (56) 

In this case, the tt pulse plays a role to accelerate the decoherence. The reason is clear from the structure of f-{t). 
When the tt pulse is applied after ii — to > 7i'/(2f^p), the two successive displacement operations in the state component 

(l T) ® DHt - h)]D[-a{h - to)] - ID® D[-a{t - h)]I)[a{ti - to)]) |n)R, (57) 

act as in-phase in the sense that the displacements by the amount of a{t — ti) and —a{ti — to) do not cancel out but 
are added in the same sign, leading to large degree of entanglement between the qubit and the reservoir. On the other 
hand, when no tt pulse is applied, the state evolves as 

(l T) ® DHt - ti)]D[a{ti - to)] + U) ® D[-a{t - ti)]b[^a{ti - to)]) \n)R, (58) 

where the two displacement operations act as out-of-phase, i.e., a(t~ ti) and a(ti — to) nearly cancel out for ti — to > 
7r/(2rip), leading to the recoherence of the qubit superposition state. 

In a realistic thermal reservoir with many modes, the above feature is more or less smeared out, but can still be 
identified at low temperatures. Roughly speaking, we had better to take the n pulse interval as 

A < 7:/{2nth), (59) 

where fith is the largest characteristic frequency of the thermalized reservoir boson spectrum /(ri)coth( ^f^). At low 
temperatures, flth is roughly fl^ in the ohmic reservoir model (Fig. or fJp in the Gaussian-ohmic reservoir model 
(Fig. 0. 



C. Inhomogeneous case 



In this subsection we take into account the inhomogeneous broadening effect in a bulk ensemble of qubits. Now we 
assume that the inhomogeneous broadening of Vj obeys the Gaussian distribution around the angular frequency oj of 
the external field 



1 



:exp 



2Sl 



(60) 
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FIG. 5: The time-resolved multiwave mixing signal intensities under the multi tt pulse irradiation (solid lines) compared with 
the ones of the single tt pulse, i.e., Af = 1 (one dotted lines) in the case of the Gaussian-ohmic resevoir model when the 
inhomogeneous broadening effect is taken into account by setting SB~5meV . The other parameters are the same as in Fig. O 



which leads us to 



^exp 



M-l 



z(z.,-c^) Am + (-1)*' ^(-1)'"A, 



dvp{i' — cij)exp 



A/-1 



i(j.-c.) Am + (-1)*' ^(-1)'"A, 



cxp 



1 AM + (-i)*^E(-ir^™pi 



(61) 



The time-resolved multiwave mixing signal intensity is finally obtained as 



it) 



P((-l)*^K(*^),t) 



-exp 



M-l 



-2r(t) - Am + (-1)*' J2 (-1)"^™ 



(62) 



The decoherence property is evaluated by measuring the echo signal. In order to see the echo signal, we consider a 
pulse sequence consisting of odd number of M tt pulses such as Aq = Ai ... = Am-2 = A, and the last interval 
Am-1 {— tu — tM-i) being varied. The echo signal is then expected around t = tM + Aa/_i as seen from Eq. I|62|l . 

Figs. 0]and[Slshows how the results of Figs. |2and|31are modified, respectively, when the inhomogeneous broadening 
effect is taken into account. We assume (5B=5meV. We see essentially the same behaviors as in Figs. |21 and |31 except 
that the echo signals are more sharpened by the fast dephasing due to the inhomogeneity. The one dotted lines 
correspond to the case of the single n pulse (M = 1), i.e., an ordinary four wave mixing scheme for photon echo 
experiment. 



IV. STATE EVOLUTION UNDER WEAK PULSES 



A. Formulation 



In an actual optical experiment, it is often the case that the maximum laser power falls short of the ideal n pulse 
and then the decoherence suppression does not work perfectly. In the qubit array system, only some fraction of qubits 
can be suppressed. In fact, the intensity of the multiwave mixing signal observed in the direction (— 1)*^K(*^^ will be 
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reduced by a factor sin^^ • • • sin''^. In other words at least this fraction of the qubits can feel the ideal tt pulses, and 
their decoherence will be suppressed. From this point of view, we consider how one can demonstrate experimentally 
the principle of the BB control. 

In the case of weak pulses, the analysis is no more straightforward. We have the complicated diffraction patterns 
in various directions. For simplicity, let us consider the three-pulse case consisting the 0th exciting pulse, the 1st and 
2nd controlling pulses. The F^^^ and G^^^ operators are given by 



(2)t 



-j<4<"' . ^0 ^2 A(2) rS(l) rS(0) i\-a^^'-°^ -a^^°>] 

= e sm— cos— cos— I)^_£)^_£)^_e ' + ' 
2 2 2 -' -' -' 

+ e-^f'cos|sin|cos|£)fJ^fJ^Se^I-^^"'+«'^'''l 
^rcoslcoslsinll^fjz^j^^i)^^'!"?"'-^^"^ 
+ e'[-*r+*^"-*f\in^sin^sin^^fJ^(^)^fJe^["^^"'+"'^"'l 



+ e" 



(63) 



and 



A(2) "U "i 

G„- = COS — cos — cos 
^ 2 2 2 



^1 jj{2} ^(1) ^(0)^,[a(."°'+c«(i°)] 



J- J- J- 



- e'^'^f -"^^"1 sin— sin— cos— DY_! D)"' e 

2 2 2 

- e^[^f'-^ricos^sin^sin^Df)Z)(i)A^°)e^ 

2 2 2 

ir0(.°)_<A<2)i . 9o 9l . r^y^i -r^y±, r^y^, 

- e ^3 'sm—cos—sm—D]_D]_/D]_/e 

2 2 2 



J210)_ 



where 



and 



ar).i:im[aP)(aW±ar)1, 



a 



(10) 



^lm[a«ari 



(64) 
(65) 

(66) 



The polarization is then given by 



P(q,f) = ^e^[i-''^+(''^-")*lTrF 



y-gi[(q-ko).r, + (.,-c.)(t-t.+Ai+Ao)]g-r___(t)lgij^^ ^Qg2^^Qg2^ 

2 2 2 

J 

^gi[(q-k,)T,+(.,-a,)(t-t.+Ai)]e-r„_o(t)i /'e^7^-o(t)cos2^ -e-^^--°Wsin2^) sin^icos^^ 



+ ^Qi[il-^2)-rj + {i^-i-u^)it-t2)]^-r.oo{t) 



1 

X- 

2 



e'^iooWcos^^ - e-^^-ooWsin^^ cos^^ + ( e^^-ooWgin^^ - e-^^-ooWcos^^ sin^^ 



2 2 



sin^o 



_ '^Qi{[<l-{-2^i-'i^o)]-rj + {''i-u,){t-t2+A^-Ao)}^-r._ + it)}_^^^Q^^^^2 
j 

_ y-ei{[q-(2k.-ko)].r, + (.,-a,)(t-t.-Ai-Ao)}e-r_++wlgjj^^ ^Qg2^^^g2^ 

^ 2 2 2 

i 

+ ^ei{[q-(2k.-k,)i.r,+(.,-a,)(t-t.-AO}e-r_+o(t)l J^_ei7-+oWcos2^ +e-^^-+''Wsin2^^ sin^isin^^ 
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gi{[q-(k2-ki+ko)]Tj+(l/j-Lj)(t-t2+Ao)} 



E 



cos7_o-(i)) — sin0osin0isin02 
g.{[q-(k2+ki-ko)].r,+(.,-^)(*-t2-Ao)} ( e-^-«+(*)cos7-o+(t)) isin0osin0isin02 



*{[q-(2k2-2ki+ko)]-r,+(i.,-u.)(t-t2-Ai+Ao)} -r_ + _(t)i • g, jjj2^gjj^2^^ 

^ 9 " 2 2 ■ 



In the above equation, we have introduced the following quantities 

rc2CiCo (t) ~ 2 



dl]/(f])Z2i£l£2|Mcoth 



fc2CiC 



C2a 



(2) 



f72 



2k^T ) ' 



(c,e{+,-,0}) 



a(2)(r!) = 


gif2t2 








— e ^ 


a(°)(r!) = 







7— oW 

7^oo(^) 
7I00W 



= 4 



= 4 



= 4 



dr2 



dr2 



m 
m 



Im 



Im 



{[a(2)(f^) + a(i)(r!)]a(°)(f^)*}, 
{[a(2)(f])_a(i)(f7)]a(")(r!)*}, 



and 



7„o-W = 7~o+W =4 



{a(2)(r!)[a(i)(f])-a(°)(r!r} 



(67) 

(68) 

(69) 

(70) 
(71) 
(72) 

(73) 
(74) 
(75) 
(76) 

(77) 



Here we are interested in the two signals observed in the directions of 2k2 — ko and 2k2 — 2ki + ko. The former 
corresponds to the photon echo signal in the four wave mixing scheme which is not affected by the first pulse (ki , 9i ) 
at ti . On the other hand, the latter corresponds to the photon echo signal in the six wave mixing scheme, and contains 
the effect caused by the (ki,0i) pulses, respectively. The signal intensities are given by 



1 



/2k2-ko(t) = 7sin"0oCOs4T!-sin^^exp[-2r_++(t) - (t - - Ai - Ao)"<5^], 



and 



1 



/2k2-2ki+ko(^) = -sin^(?osin^ysin4^exp[~2r_+_(t) - (t - ^2 - Ai + Ao)"<5|], 



(78) 
(79) 



where the inhomogeneous broadening distribution of Eq. H6U|) has been assumed. 

In the following we assume that 6q = 9i = 62 = n /"i. This pulse area is recently achievable in the laboratory for 
common semiconductor quantum well and dot systems, such as GaSe and GaAs. Note that the signal intensities 
observed at 2k2 — ko and 2k2 — 2ki +ko result from the qubits that effectively evolve under the (ko, 7r/2 ; ki, tt) and 
(ko, ■7r/2 ; ki, tt ; k2, tt) pulse sequence. The two signal intensities have the same weight factors. So we can directly 
compare the signal intensities observed at the two directions /2k2-ko(*) and /2k2-2ki+ko(*) with each other. The 
difference between them is due to the effect by the first tt pulse, which is what should be confirmed experimentally. 
One suitable way to see this difference is to sweep ti by fixing t2 such that the signal intensity /2k2-ko(0 remains 
constant while /2k2-2ki+ko (i) varies. 
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FIG. 6: The intensities of the time-resolved six wave mixing signals in the case of homogeneous qubits and the ohmic reservoir 
model with a = 0.1, Qc = 4meV,and T = 100 K. The thick solid line corresponds to the signal intensity /2k2-2ki+ko (i)) while 
the thick dashed line to /2k2-ko(i)- The one dotted line represents the free induction decay signal. The time of the second 
pulse is added at t2 — 0.2 ps. The time of the first pulse ti — 0.02, 0.04, and 0.18 ps is indicated by the vertical dashed line. 



B. Homogeneous case 



Let us start with the case of homogeneous qubits {Sb — 0). Fig. shows the signals /2k2-2ki+ko (0 (the sohd 
hne) and /2k2~ko(^) (the dashed line) for various values of ti setting t2 — 0.2 ps. We have set to = hereafter. As a 
reference, the free induction decay signal /kf, (t) is also shown by the one dotted line. The reservoir model is assumed 
to be the ohmic model with a = 0.1, ilc — 4meV,and T — lOOK. The second pulse is added at t2 = 0.2ps. The 
time of the first pulse ti = 0.02, 0.04, and 0.18 ps is indicated by the vertical dashed line. The effect of decoherence 
suppression by the (ki,6'i) pulse is seen for ti = 0.02, 0.04, and 0.12ps as the difference between the solid and 
dashed lines. 

When the signal /2k2-2ki+ko(i) is greater than /2k2-ko('^) for almost all regions of t > t2, (the cases of ti — 
0.1, 0.2, and 0.6ps in Fig. 0, one may use the time-integrated signal intensity, i.e., the total energy of the emitted 
radiation 

r^' = / dti^it) (80) 

to measure the effect of decoherence suppression. In fact, the measurement of the time-integrated signal intensity is 
much easier than that of the time resolved signal intensity Iq{t). Fig. |7| shows the time-integrated signal intensity 



t, = 0.2 (ps) = 0.4 (ps) = 0.6 (ps) 




f, (ps) f, (ps) f, (ps) 



FIG. 7: The time-integrated signal intensities as a function of the first pulse ti with t2 fixed. The homogeneous qubits and 
the ohmic reservoir model are assumed. The details are explained in the text. 



as a function of <i, fixing the second pulse time t2- The upper three viewgraphs correspond to the case of 

the ohmic reservoir with a — 0.1 and Sic — 4meV. ^2=0. 2, 0.4, 0.6 ps, from the left. Each viewgraph contains the 
cases of three kinds of temperature T =10, 50, and 100 K. The lower three viewgraphs are the same as the upper but 
the higher cutoff frequency flc = 8 meV. Since the signal at the direction 2k2 — kg is not affected by the first pulse, 
the time-integrated signal intensity as a function of ti, -f2k2-k;o(^i)' results in a flat horizontal line. On the other 
hand, -f2k2-2ki+ko (^i) shows peak structure at certain time(s) of ti. For the lower cutoff frequency fie and higher 
temperature, ^2k2-2ki+ko (^i) becomes larger than -^2ic2~ko(^i) ^^'^ peaked around t — t2/2 or an earlier time. 

For the higher cutoflF frequency fie and lower temperature (the case of 10 K in Figs. [7^ and f), on the other hand, 
there appears a time region of ti where ^2k2-2ki+ko (^i) becomes smaller than -f2icQ-ko (^i)' ^-^^ ^^'^^ pulse (ki,6'i) 
fails in suppressing the decoherence and rather accelerates it. The origin of this phenomenon is the in-phase coupling 
between a qubit and the reservoir bosons caused by the first (intermediate) pulse (ki,0i), as explained in section 
1111 B\ In order to suppress the decoherence effectively, the pulse interval should be A < 7r/(2rith), where Slth is the 
characteristic frequency at which the thermalized boson factor defined by 

^in,T)^i{n)cothi^), (81) 

is peaked. Note that the decoherence exponent in Eq. (|39() can be expressed as 

r{t)= / dnan,tun,T), (82) 

using j]{n,T) and the time dependent factor defined by 

an,t)^2l^. (83) 



The condition A < 7r/(2r2th) is not satisfied in the case of 10 K in Figs. [7|3 and f. In fact, for the higher cutoff 
ferequency Qc and lower temperature, Qth becomes higher as seen in Fig. |H1 
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2.5 
> 2.0 




n (meV) 

FIG. 8: The thermalized boson factor r;(f2, T) in the case of the ohmic reservoir model with f2c=8meV and a = 0.1. The five 
curves correspond to the temperatures, from the bottom, IK, lOK, 20K, 5QK, and lOOK. 

It is commonly encountered in semiconductor exciton systems that the decoherence is accelerated when an additional 
pulse is irradiated. This is usually due to that an additional pulse increases the number of excitons, which in 
turn enhances the exciton-exciton interaction, and increases the number of decoherence channels. In the present 
model, however, the interaction between the qubits themselves is not taken into account. Therefore the decoherence 
acceleration effect purely comes from the coherent interaction between a qubit and the reservoir bosons. More 
remarkably in the case of lOK in Figs. 0: and f, ^2k2-2ki+ko (^i) exceeds ^2k2-ko(^i) ^-gain, and forms a peak for 
larger t^. In this region, the out-of-phase coupling between a qubit and the reservoir bosons dominates, and the first 
pulse (ki, 6'i) succeeds in suppressing the decoherence. This kind of cross over from the decohernce acceleration to the 
suppression is a peculiar aspect of the present model, which is expected at low temperatures, and must be interesting 
in its own light to be investigated experimentally. 

C. Inhomogeneous case 

Now let us consider the inhomogeneous broadening effect in a bulk ensemble of qubits. Fig. El shows how the 
time-integrated signal intensities change when the inhomogeneous broadening effect is introduced into the ohmic 
reservoir model. From the upper viewgraphs the inhomogeneous broadening effect is taken as (5b=0, 2, and 5meV. 
The reservoir parameters are a = 0.1 and ^lc = 8meV. The time of the second pulse ^2 is 0.2 ps and 0.6 ps for the left 
and right viewgraphs, respectively. It should be noted that the inhomogeneously broadened signals /2k2-2ki+ko (t) and 
^2k2-ko(0 appear as the echos around the times ^ = ^2 + ^1 — ^0 =2(^2 — ^1) and i = ^2 + ^i + Aq = 2^2, respectively. 
While the temporal position of the echo signal hk^-kgit) does not depend on ti, the echo signal /2k2-2ki+ko(i) appears 
earlier as ti gets longer delay. In particular, when the time ti exceeds ^2/2 such that Ai < Aq, the echo signal reduces 
exponentially. As seen, it is still possible to observe the decoherence suppression effect, i.e. that ^2k2-2ki+ko (^1) 
becomes larger than -?^2k2-ko(^i) even in the presence of the inhomogeneous broadening effect. On the other hand, the 
decoherence acceleration effect seen in the right top viewgraph is rapidly smeared out as the amount of inhomogeneous 
broadening increases. 

The decoherence acceleration effect can be observed more clearly in the Gaussian-ohmic reservoir model, which has 
a sharper spectral density around a characteristic frequency fip than that in the ohmic reservoir model. Fig. llUl shows 
the thermalized boson factor ry(il,r). Fig. 1111 shows the time integrated signal intensities for the Gaussian-ohmic 
reservoir model with the parameters ftp =13meV, 7p =4meV, ap = 0.05, flc =8nieV, and a = 0.1. From the upper 
viewgraphs the inhomogeneous broadening effect is taken as 6b =0, 2, and 5meV. The time ^2 is 0.2 ps and 0.6 ps for 
the left and right viewgraphs, respectively. Among the six viewgraphs, we would like to pay attension to Figs. Illb . 
e, and f. In this case, -^2k*,-2ki+ko(^i) exhibits the decoherence acceleration effect first, the decoherence suppression 
effect then, as a function of ti. Such a cross over is a very characteristic evidence of the coherent qubit-reservoir 
interaction in the non-Markovian region. More importantly, this behavior can be seen even under the realistic amount 
of the inhomogeneous broadening effect such as Sb =5meV. This may appeal to experimental investigation on the 
coherent qubit-reservoir interaction. 
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FIG. 9: The time-integrated signal intensities as a function of the first pulse ti with t2 fixed. The ohmic reservoir model is 
assumed with the parameters a — 0.1 and Qc =8meV. The top figures correspond to the homogeneous case, while the middle 
and bottom ones to the inhomogeneous case with 5B~2meV and 5B~5meV , respectively. 
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FIG. 10: The thermalized boson factor 7?(n, T) in the case of the Gaussian-ohmic reservoir model with Qp =13 meV, 7p =4meV, 
ftp — 0.05, flc =8meV, and a = 0.1. The five curves correspond to the temperatures, from the bottom, IK, lOK, 20K, 50K, 
and lOOK. 



17 



t= 0.2 (ps) 



t= 0.6 (ps) 



w 

• 

c 

3 

«, 

>< 

'5) 
c 

0) 
+.« 

c 

« 
c 

D) 

S) 

T5 
0) 



0) 



0) 

E 



0.10 
0.05 


0.04 
0.03 
0.02 
0.01 
0.0^ 
0.015 

0.010 

0.005 

0.01 




.00 0.05 0.10 0.15 0.20 
t, (ps) 



8,=2meV 2k^-2K+k„,^ 



00 0.05 0.10 0.15 0.20 
t, (ps) 




r.00 0.05 0.10 0.15 0.20 
t, (ps) 



0.03 
0.02 
0.01 
0.0| 

0.004 

0.002 

O.OOg 
0.004 

0.003 
0.002 
0.001 
0.001 



8^=0 meV 2k -2k +k 




1.1 0.2 0.3 0.4 0.5 0.6 
^(ps) 



e 

8„=2meV ;^ 

/ \ 

/.X I 


t,-2K+k„ 

2 10 

2k-k„ 

2 







1.1 0.2 0.3 0.4 0.5 0.6 
^(PS) 




'8.0 0.1 0.2 0.3 0.4 0.5 0.6 

t, (ps) 



FIG. 11: The time-integrated signal intensities as a function of the first pulse ti with t2 fixed in the case of the Gaussian-ohmic 
reservoir model with the parameters f2p =13meV, 7p =4meV, ap = 0.05, fie =8meV, and a = 0.1. The top figures correspond 
to the homogeneous case, while the middle and bottom ones to the inhomogeneous case with (5B=2meV and 5B=5meV, 
respectively. 



V. CONCLUDING REMARK 



We have developed a theory to analyze the decoherence in the qubit array system with the photon eeho signals in 
the multiwave mixing configuration. We have presented how the decoherence suppression effect by the BB control 
with the TT pulses can be demonstrated in laboratory by using a bulk ensemble of excitons and optical pulses whose 
pulse area is even smaller than tt. The key is to analyze the time-integated multiwave mixing signals diffracted into 
certain phase matching directions from a bulk ensemble. We have given numerical examples in the six wave mixing 
configuration (three pulses), where the signals in two directions 2k2 — ko and 2k2 — 2ki + ko are compared with each 
other by sweeping the pulse intervals. Depending on the pulse interval conditions, both the decoherence suppression 
and acceleration effects take place as pointed out in earlier works. The cross over from one to the other is a clear 
evidence of the coherent qubit-reservoir interaction. We have shown that this cross over may be observed even under 
realistic inhomogeneous broadening. This encourages experimental investigations in qubit systems in solid state. 

To understand this cross over, we have introduced the notions of in-phase and out-of-phase couplings between a 
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qubit and the reservoir bosons. The decoherence process in the present model is mathematically described by the state 
evolution caused by the displacement operations on the reservoir boson modes, depending on the qubit states of 
components. The in-phase coupling means that the successive displacement operations between the optical pulses act 
additively, enhancing the entanglement between a qubit and the reservoir bosons. The out-of-phase coupling means 
that the successive displacement operations act as canceling out. 

Which type of coupling dominates is determined by the relation between the pulse interval and the reservoir boson 
spectrum. Therefore by analyzing such behaviors systematically one may identify the reservoir characteristics. In 
particular, when analysis is made in the weak pulse configuration where we can observe many diffraction signals 
from various multiwavc mixing channels simultaneously, we may increase the precision of parameter estimation by 
fitting as many diffracted signals as possible with a single theoretical model. Thus our theory will be useful to 
investigate the reservoir characteristics, which is an important first step to proceed quantum information processing 
with semiconductor excitons. 
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